Vortices in Bose-Einstein condensates with anharmonic confinement 
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We examine an effectively repulsive Bose-Einstein condensate of atoms, that rotates in a quadratic- 
plus-quartic trapping potential. We investigate the phase diagram of the system as a function of 
the angular frequency of rotation and of the coupling constant, demonstrating that there are phase 
transitions between multiply- and singly-quantized vortex states. The derived phase diagram is 
shown to be universal and exact in the limits of small anharmonicity and weak coupling constant. 
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The behavior of Bose-Einstein condensates of alkali- 
metal atoms under rotation has become a very inter- 
esting topic in recent years, since these systems are su- 
perfluid and the only way to rotate them is by creating 
quantized vortex states. The experimental observation 
of vortex states in a two-component system was reported 
by Matthews et al. [1], while Madison et al. [2] observed 
vortex states in a stirred one-component Bose-Einstein 
condensate. Many more experimental studies on rotat- 
ing condensates have been performed recently [3-5]. 

If, in a harmonic trapping potential, the frequency of 
rotation becomes equal to the trap frequency, the cen- 
trifugal force exactly cancels the restoring force and the 
atoms fly apart. In other words, the trap frequency 
sets an upper limit for the frequency of rotation of the 
cloud. To overcome this difficulty and also to examine 
the various phases which have been predicted theoreti- 
cally [6-13], Bretin et al. [14] recently examined exper- 
imentally a rotating Bose-Einstein condensate in an an- 
harmonic potential having the form 
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Here, p is the cylindrical polar coordinate, M is the 
atomic mass, ui is the frequency of the harmonic trapping 
potential, ao = {Ti/Muj) 1 / 2 is the oscillator length, and 
A is a dimensionless constant which in the experiment of 
Ref. [14] was very small, on the order of 10~ 3 . 

In what follows we investigate the behavior of a gas 
that is trapped in a potential of the form of V and is 
under rotation. In reality, the gas is also trapped in the 
z direction. However since it rotates around this axis, we 
neglect the trapping in this direction and assume that 
the gas has a constant density per unit length a = N/Z, 
where N is the number of atoms and Z is the width of the 
condensate. In addition we assume that the degrees of 
freedom along the z axis are frozen so that our problem 
becomes two-dimensional. In this case and for a purely 
harmonic potential, A = 0, the single-particle energy lev- 
els are given by 



where n r is the radial quantum number, and m is the 
quantum number corresponding to the angular momen- 
tum. Since e grows linearly with the angular momen- 
tum, any superposition of states with different values of 
m (and total angular momentum mo) will have the same 
energy as the pure state with m — mo- Interactions split 
this degeneracy. As shown in Refs. [15,16], repulsive in- 
teractions always favor a superposition of states. The ar- 
gument is very much like that in the tight-binding model, 
where derealization of the wavefunction lowers the en- 
ergy. In practice, this implies that a singly-quantized 
vortex state always has a lower energy than a multiply- 
quantized vortex in a harmonic potential. 

This picture changes drastically when A > 0. The 
single-particle energy levels grow faster than linearly with 
m in this case. As a result, for sufficiently weak interac- 
tions, a multiply-quantized vortex state has lower energy 
than a singly-quantized vortex with the same circulation. 
As the coupling constant increases, however, it eventually 
becomes energetically favorable for the system to spread 
its angular momentum to additional states, and multiply- 
quantized vortex states break into singly-quantized vor- 
tices. We demonstrate this phase transition explicitly 
below. 

In the following, we assume that the interaction is de- 
scribed by a contact potential of the form 
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e = (2n r + \m\ + l)hui, 



(2) 



where Uq — ATrh 2 a/M is the strength of the effective two- 
body interaction with a equal to the scattering length for 
atom-atom collisions. We assume that the interaction is 
repulsive, a > 0. Attractive interactions in anharmonic 
potentials have been studied in Ref. [11]. 

We can estimate the interaction strength required at 
the transition for a given value of the anharmonicity co- 
efficient, A, by noting that the characteristic energy scale 
associated with the anharmonicity must be comparable 
to the interaction energy. (This argument, valid when 
the number of vortices is of order 1, will be generalized 
below.) The first scale is of order Xhuj, and the second 
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is of order anllo, where n is the characteristic atom den- 
sity, and a is a dimensionless constant of order 1CP 1 . 
(See, e.g., Fig. 1 of Ref. [16].) If Xhuo is to be w anU , 
the effective coupling constant, aa, must be w A/a. For 
A <~ 10~ 3 , we conclude that aa must be less than 1. 

Therefore, provided that A and aa are smaller than 
unity, we consider the cigenstates of the harmonic po- 
tential with zero radial excitations and angular momen- 
tum mh (to is positive for simplicity) as our unperturbed 
states, 
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where <j) is the angle in cylindrical polar coordinates. The 
expectation value of K E + V in the state $ m , where KE 
is the kinetic energy, —h 2 V 2 /2M, is 



($ m \KE + V\$ 

raj 



huj[m + ^(m + l)(m + 2)] (5) 



measured with respect to the zero-point energy, hui. (All 
energies will henceforth be measured with respect to hu>.) 
Equation (5) shows explicitly the faster-than-linear (i.e., 
quadratic) dependence of the energy on to. 

Let us now incorporate the interactions. We use per- 
turbation theory to determine the interaction energy per 
particle in the state $ m , 
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Combining Eqs. (5) and (6) we obtain the total energy 
per particle in <!> m , 
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A, 1W „ N (2m)! 
+ -(m + l)(m + 2) + aa ^_, 
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In a rotating frame of reference with some angular veloc- 
ity Q, the energy is 



N hto 



m(l--) + -(m + l)(m + 2) + < ra^ i5 5. 

(8) 



The above expression gives the energy per particle in the 
rotating frame in terms of the parameters Q,/uj, A, and 
era, in a given state $ m . 

From Eq. (8) we see that E[ becomes lower than E' 
for a critical frequency of rotation which is 



flx/w = 1 + 2A - aa/2. 



(9) 



This is actually the critical angular frequency for induc- 
ing rotation in the cloud. For A = 0, this result is in 
agreement with Refs. [15,16]. As fl increases further, 
states with larger m become energetically favorable. The 
critical value of fl m for a multiply-quantized vortex state 
to form with mh units of angular momentum is given by 
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FIG. 1. The angular momentum, m, as function of the ro- 
tational frequency, Q,/uj, as given by Eq. (10) with A = 0.05 
and aa — 0.4. As fl/iv increases, multiply-quantized vortices 
form at the critical frequencies indicated in the graph. 



flm/ui = 1 + A(m + 1) — aa 



(2m- 2)! 



2 2m " 1 (m- l)!mf 



(10) 



A similar result was derived using a variational approach 
in Ref. [11] for an attractive Bose gas trapped in all three 
directions. According to our model and for sufficiently 
weak interactions, the angular momentum mh of the sys- 
tem is quantized and increases stepwise at the frequen- 
cies given by Eq. (10), as shown in Fig. 1 for A = 0.05 and 
aa = 0.4. For m » 1, fl m /Ld w 1 + Am. Figure 2 shows 
the phase diagram of multiply-quantized vortex states in 
the n/cu - aa plane, as given by Eq. (10), for A = 0.05. 

We argued earlier that, as the strength of the inter- 
action increases, a multiply-quantized vortex state will 
split into singly-quantized vortices. To attack this prob- 
lem more generally, one can expand the order parameter 
[15,16], 



(11) 



Here the coefficients c„ 
eters subject to 



are (complex) variational param- 



TO C™ = 



(12) 



The first condition in simply the normalization, while the 
second implies that the state ^ is assumed to have some 
angular momentum per particle, lh, where / is not nec- 
essarily quantized. After minimization of the energy in 
the rotating frame with respect to the variational param- 
eters c m and with respect to I, one can determine \& and 
I = l(Sl). 
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The phase transition we consider first appears in the 
case with m = 2, and this is the example we consider 
below. However our method can be generalized to higher 
values of m. Let us therefore consider the trial wavefunc- 
tion [15,16] 



Co$0 + C 2 $2 + C 4 $4- 



(13) 



The state VP reduces to $2 when \c2\ = 1, Co = C4 = (i.e., 
it is a doubly-quantized vortex state), and it describes 
two singly-quantized vortices when c c 4 7^ [17]. As we 
show below, there is a critical value of aa above which Co 
and C4 become nonzero. It is important to note that, as a 
consequence of the spherical symmetry of the interaction, 
the phase boundary between single and multiple quanti- 
zation is given exactly by the trial order parameter of 
Eq. (13) so long as the states with nonzero radial nodes, 
n r 7^ 0, can be neglected. Equation (13) can be extended 
to larger values of to by considering the local stability of 
the multiply-quantized state, $ m , with respect to the ad- 
mixture of states with any angular momentum mi and 
to 2 = 2m — mi. For small to, the leading instability is 
given by mi = 0. These results will be addressed in a 
future publication. 

Calculating the energy per particle in the state W in 
the rotating frame of reference we get, 



1_£? 



= 1(1 - Sl/w) + A(|c | 2 + 6|c 2 | 2 + 15|c 4 | 2 ) 
-aa(|co| 4 + | C o| 2 | C 2| 2 + ^| C2 | 4 + J| C o| 2 |c4| 2 
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+ 128 N + 16 |C2 ' N " X |ColN N) - (W) 

Using Eqs. (12), we can eliminate |co| and | C4 1 in Eq. (14). 
For a fixed A, we then vary C2 and I to minimize E' for 
various values of Ct/u) and a a and thus investigate the 
phase diagram in the fi/w - aa plane. The result for 
A = 0.05 is shown as the dashed line in Fig. 2 in the m = 2 
phase. Above this line there are two singly-quantized 
vortices with |c2 1 < 1- Below this line there is a doubly- 
quantized vortex state, where |c 2 | = 1. The actual phase 
boundary is expected to lie lower than the dashed line be- 
cause of the variational nature of the present calculation 
(i.e., the neglect of states with n r 7^ 0). Extending the 
same approach to higher values of m, we derive the cor- 
responding phase boundaries for m = 3, 4 and 5, which 
are shown as dashed lines in Fig. 2. The overall phase 
diagram of Fig. 2 resembles that obtained numerically by 
Lundh [6] for traps of some power law. 

It is important to note that the transition seen here is 
continuous and of second order. Thus, it is sufficient to 
consider the local stability of the states $ m in determin- 
ing the phase boundary. As aa increases, the separation 
between the two vortices increases continuously in agree- 
ment with numerical simulations [6]. Note also that the 
form of the phase diagram is unaltered by changes in the 
strength of the anharmonicity since (1 — Q/uj) and aa at 
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FIG. 2. The phase diagram of vortex states in a 
quadratic-plus-quartic potential in the Q/to - aa plane for 
A = 0.05. The solid lines are the phase boundaries between 
multiply-quantized vortex states with in = 0, 1, 2, 3, 4, 5, and 
6, as given by Eq. (10). The dashed lines indicate the bound- 
ary between the phase of singly-quantized vortices, and the 
phase of multiply-quantized vortex states. 



the phase boundary simply scale linearly with A. The 
phase diagram shown in Fig. 2 is thus universal and ex- 
act in the limits of small anharmonicity and weak inter- 
actions. 

Another important point is that in the phases of sin- 
gle quantization, the function I = 1(0) is not necessarily 
quantized. Therefore, the graph of Fig. 1 for higher val- 
ues of aa will be partly continuous and partly discontinu- 
ous. This is similar to the case of harmonic confincmemt 
[15] except for the fact that, in the present problem I 
does not diverge when O — > u) because of the presence of 
the quartic term in the trapping potential. 

To understand better the mechanism which gives rise 
to the phase transition between single and multiple quan- 
tization, we return to the trial order parameter of Eq. (13) 
and consider the extra constraint 1 = 2. Eliminating |co| 
and \ci\ in the energy, we find in the laboratory frame 
that 



1 E 



2 + 2A(4- |c 2 | 2 ) 
(64V6- 106)|c 2 p 



+ (64V6- 109)|c 2 | 4 ]. (15) 



The right side of the above equation must be minimized 
with respect to |c 2 |. However, it is clear that the en- 
ergy that results from the anharmonicity is minimized 
for I C2 1 = 1) while the interaction energy is minimized 
for |c 2 | 2 = (64V6 - 106)/2(64V6 - 109) w 0.53, which 
is less than 1. It is essentially the balance between these 
two terms that gives rise to the phase transition between 
multiple and single quantization. Some trivial algebra 
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shows that the critical coupling constant above which 
the vortices separate is given by 



(aa) c = 



64A 



4\/6 - 7 



22.88A. 



(16) 



The dashed curve in Fig. 2 (for m = 2) saturates for 
larger f2 to the value of aa given by Eq. (16). (In fact, it 
has essentially reached this point at the highest value of 
aa shown.) Since the solid line separating the m = 2 and 
m = 3 phases in Fig. 2 cuts this dashed line at some aa < 
(aa) c , Eq. (16) provides an upper limit for aa for which 
there is a doubly-quantized vortex state. In Ref. [6] the 
function (aa) c — (aa) c (X) was calculated numerically (in 
Fig. 4 of [6]), and (aa) c of Eq. (16) indeed lies above this 
curve. The discrepancy is less than 10 percent for the 
smallest A shown, A = 0.01, and grows as A increases. 

We note that Ref. [10] has predicted the existence of a 
triple point in the Cl/u) - aa plane, at which a multiply- 
quantized vortex state, a vortex lattice, and a vortex lat- 
tice with a hole coexist. The present study is far from 
this point. However, such a triple point would be in- 
dicated in our approach when the leading instability is 
due to the admixture of states with mi ^ 0, since $ is 
the only component in the wavefunction that does not 
vanish at p = 0. Since the effective trapping potential, 
V — Mtt 2 p 2 /2, has the familiar "Mexican hat" shape for 
Sl/w > 1, the contribution of $0 to \I> will decrease with 
increasing Cl/co. 

Turning to the experiment of Ref. [14], for A = 10~ 3 
Eq. (16) implies a value of (aa) c 0.02, which is very 
small. Actually this experiment is closer to the Thomas- 
Fermi limit of strong interactions [18]. For a number of 
atoms N = 3x 10 5 , u = 2ttx65.6 Hz, lu z = 2ttx 11.0 Hz, 
and a = 53 A, we find a transverse radius of the cloud 
of w 6.6 /im, and a width in the z direction of ss 39.4 
/im. With these values aa turns out to be w 20, which 
is much greater than (aa) c . This fact explains why no 
multiply-quantized vortices were observed. On the other 
hand, tens of vortices were created in this experiment, 
and (aa) c increases with the number of vortices [10]. We 
can obtain a more general estimate for (aa) c for the m-th 
multiply-quantized vortex state as 



(aa) c — A- 



" 1 m!(m + 2)! 
2m! : 



(17) 



by equating the energy due to the quartic term, Eq. (5), 
with the interaction energy, Eq. (6). Equation (17) is in- 
deed an increasing function of m. In Eq. (17) A(m + 
l)(m + 2)/2 must be smaller than m, and m should 
thus not exceed 2/A. Self-consistency requires that (aa) c 
should not exceed unity. To consider a specific exam- 
ple with A = 10~ 3 for a multiply-quantized vortex state 
with m = 15, Eq. (17) implies that (aa) c ~ 0.94, while 
Eq. (10) gives a rotation frequency ^15 w l.Olw. 

In summary, we have examined a weakly-repulsive 
Bose-Einstein condensate of atoms that is confined in 
a quadratic-plus-quartic trapping potential and is under 



rotation. We have investigated the two phases that the 
system exhibits - multiply- and singly-quantized vortices 
- as functions of the frequency of rotation, the coupling 
constant and the anharmonicity. We have also demon- 
strated that, by varying the coupling constant, the cloud 
undergoes a phase transition from multiply- to singly- 
quantized vortex states. 
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